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Abstract 

The present study was undertaken to investigate growth kinetics and time-dependent change in global expression of 
Escherichia coli 0157:H7 Sakai upon an abrupt downshift in water activity (a„). Based on viable count data, shifting £ coli 
from a„ 0.993 to a„ 0.985 or less caused an apparent loss, then recovery, of culturability. Exponential growth then resumed 
at a rate characteristic for the a^,„ Imposed. To understand the responses of this pathogen to abrupt osmotic stress, we 
employed an Integrated genomic and proteomic approach to characterize its cellular response during exposure to a rapid 
downshift but still within the growth range from 0.993 to a^„ 0.967. Of particular Interest, genes and proteins with cell 
envelope-related functions were Induced during the initial loss and subsequent recovery of culturability. This Implies that 
cells undergo remodeling of their envelope composition, enabling them to adapt to osmotic stress. Growth at low a^,„, 
however, involved up-regulating additional genes and proteins, which are involved In the biosynthesis of specific amino 
acids, and carbohydrate catabolism and energy generation. This suggests their Important role In facilitating growth under 
such stress. Finally, we highlighted the ability of £ coli to activate multiple stress responses by transiently Inducing the RpoE 
and RpoH regulons to control protein misfolding, while simultaneously activating the master stress regulator RpoS to 
mediate long-term adaptation to hyperosmolality. This investigation extends our understanding of the potential 
mechanisms used by pathogenic £ coli to adapt, survive and grow under osmotic stress, which could potentially be 
exploited to aid the selection and/or development of novel strategies to inactivate this pathogen. 
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Introduction 

Enterohemorrhagic Escherichia coli (EHEC) 0157:H7 have 
emerged as important food-borne pathogens of considerable 
public health concern. They can cause a range of human illnesses 
including diarrhea, hemorrhagic colitis, and the life-threatening 
hemolytic uremic syndrome [1,2]. The majority of reported 
outbreaks and sporadic cases of 0157:H7 infection appear to be 
attributed to the consumption of foods of bovine origin, although 
cases involving dairy products, water, vegetables and fruit products 
have also been reported. Numerous studies have also identified 
ruminant animals, especially cattle, as the major reservoir of £. coli 
0157:H7, which is usually found in the faeces and rumen, on the 
hide and derived carcass surfaces [1-3]. 

Increasing the osmotic pressure (lowering a„) is one of the most 
widely used methods in food preservation to control the growth of 
bacteria, including E. coli. Reduction in the external a„ typically 
results in a rapid loss of the cytoplasmic volume in a process called 
plasmolysis, and causes reduced respiration and growth arrest, 
whereas both intracellular ATP and cytoplasmic pH have been 
reported to increase [4,5]. To adapt to hyperosmotic stress, 
bacteria employ adaptive mechanisms referred to generally as 



osmoregulatory systems. A major role of these systems is to 
maintain the proper intracellular osmotic pressure within tolerable 
limits. This generally involves accumulation of charged solutes 
(e.g. potassium ions (K"*^ and glutamate), followed by accumulation 
of compatible solutes either through de now biosynthesis (e.g. 
trehalose) or through uptake from the external environment (e.g. 
glycine betaine and proline) [6,7]. Furthermore, it is well 
established that bacterial cells previously exposed to osmotic stress 
(i.e. osmoticaUy-adapted cells) acquire increased resistance to other 
stresses such as high temperature and oxidative stresses [6,8]. 
Therefore, the ability of pathogenic bacteria to adapt to and 
survive under adverse conditions could increase the risk of food- 
borne illness. A detailed understanding of how E. coli 0157:H7 
adapts to hyperosmotic stress could aid in identification of 
potential targets and develop effective interventions for controlling 
or eliminating this pathogen. 

Previously, we employed both cDNA microarray and 2D-LC/ 
MS/MS analyses to elucidate the genome and proteome 
expressions of exponential phase E. coli 0157:H7 strain Sakai 
grown under steady-state conditions, relevant to low temperature 
and water activity conditions experienced during carcass chilling 
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[9]. It was found that E. cofi 0157:H7 respond to these steady-state 
conditions, including osmotic stress (25°C ii„ 0.967) by activating 
the master stress response regulator RpoS and the Res phosphor- 
elay system involved in the biosynthesis of the exopolysacharide 
colanic acid, as well as down-regulating genes and proteins 
involved in chemotaxis and motility [9]. Such findings have 
provided a basehne of knowledge of the potential molecular 
mechanisms enabling growth of this pathogen under these stress 
conditions. To gain a deeper insight into the physiology of 
exponentially growing E. coli 0157:H7 Sakai in response to 
hyperosmolality, the present study investigated the growth kinetics 
of this pathogen subjected to sudden osmotic upshift (from ii„ 
0.993 to a^ 0.980, 0.975, 0.970, 0.967 or 0.960), as well as to 
examine the time-dependent alterations in its transcriptome and 
proteome upon hyperosmotic shock from a„ 0.993 to a„ 0.967 at a 
constant temperature of 35°C. The global responses were analyzed 
by botii cDNA microarray and 2D-LC /MS/MS analyses. 
Differences in gene and protein expression patterns in E. coli 
before and after osmotic shock were analyzed through quantitative 
and comparative analysis of time-series changes in both mRNA 
and protein levels. 

Materials and Methods 

Division of experimental worl< 

Parallel transcriptomic and proteomic studies were undertaken 
in two independent laboratories, CSIRO Animal, Food and 
Health Sciences (North Rydc, Australia) and Food Safety Centre 
(Hobart, Australia), respectively. A series of experiments on the 
growth response of E. coli 0157:H7 to hyperosmotic shifts were 
also conducted in the Food Safetj- Centre (Hobart, Australia). 
Although conducted in separate laboratories, the protocol was 
identical for the hyperosmotic shock experiments with the only 
variation involving the different methods used for extraction of 
mRNA and proteins, as well as the use of different equipment for 
optical density measurement (Transcriptomic: Pharmacia Biotech 
Novaspec 11 spectrophotometer; and Proteomic: Spectronic 20, 
Bausch and Lomb) and waterbaths for incubation of cultures 
(Transcriptomic: Julabo SW23 waterbath (80 rpm); and Proteo- 
mic: Ratek SWB20D shaking waterbath (80 osciUations/min)). 

Bacterial strain and preparation of inocula 

E. coli 0157:H7 strain Sakai [10,1 1] obtained from Dr. Carlton 
Gyles (Guelph, Canada) was used in this study. The strain was 
maintained at — 80°C. Prior to any experiments, the culture was 
revived by plating onto brain-heart infusion (BHI) agar (Oxoid, 
CM225), and incubating at 37°C for 24 h. A well-isolated single 
colony was aseptically transferred into 2.^ ml of BHI broth. After 
incubation at 37°C for 20 h, the broth culture was stored at 4°C 
and used as a 'working' culture within a week. 

Preparation of low a^ broth 

A BHI broth containing a very high concentration of NaCl was 
prepared to be used in osmotic upshift experiments. Specifically, 
an appropriate amount of NaCl was added to fresh BHI broth (a„ 
0.993±0.00,3) to achieve the final a„ of 0.760 (±0.003), as 
determined from the Tables of Chirife and Resnik [12]. This 
NaCl-concentrated BHI broth was then sterilized (121°C, 20 min) 
and kept at room temperature until used. The final a„ of each 
NaCl-modified BHI broth was confirmed using an Aqualab CX-2 
dew point instrument (Decagon Devices, Inc., Pullman, USA). 



An abrupt downshift in water activity 

An appropriate volume of the working culture was diluted 1:10* 
in 25 ml of pre-warmed (35°C) BHI broth. This 'primary' culture 
was incubated at 35°C in a water bath with agitation (80 oscil- 
lations/min; Ratek Instruments, Boronia, Australia), and its 
growth monitored turbidimetricaUy at 600 nm with a Spectronic 
20 spectrophotometer (Bausch and Lomb, Rochester, USA). After 
achieving an ODroo of 0. 1 ± 0.0 1 (i.e. the mid-exponential phase of 
growth, —10^ CFU/ml), a 1:10^ dilution of the primary culture 
was prepared further in a 25 ml of pre-warmed (35°C) BHI broth 
to prepare a 'secondary' culture. This culture was then re- 
incubated at 35 °C with agitation until it reached the mid- 
exponential growth phase (ODgoo of 0.1 ±0.01). At this point, after 
removing an aliquot for subsequent analysis, an appropriate 
volume of the NaCl-concentrated BHI previously prepared (pre- 
warmed to 35°C) was added to the culture to create a series of 
hyperosmotic brotiis at a„ 0.980, 0.975, 0.970, 0.967 or 0.960. 
Each of these inoculated ajv-modified broths was re-incubated at 
35°C with shaking. 

Microbiological analysis 

E. coli growth in each a^^-modified broth was monitored by 
viable count method. At regular intervals, 0. 1 ml of the test broth 
was obtained and serially diluted in 0. 1 % peptone water (PW; 
Oxoid, LP0037) containing 0.85'/<. NaCl. Aliquots of the 
appropriate dilution were then plated onto BHI agar supplement- 
ed with 0. 1 % sodium pyruvate (Sigma, P8574) using an Autoplate 
4000 spiral plater (Spiral Biotech, Bethesda, USA). AH plates were 
then incubated at 37°C for 24 h before colonies were counted. All 
growth experiments were rephcated at least twice, and the results 
reported as the average of the replicates. 

Logio CFU/ml was plotted against time to generate growth 
curve for each osmotic treatment. AH growth curves were analyzed 
using linear regression to estimate generation time (GT) and lag 
time (LT), according to Mellefont et al. [13]. Briefly, a straight line 
was fitted to the data points (i.e. logio CFLV ml) that best represent 
exponential phase of growth. To ensure the accuracy of this 
estimation, at least four time points over this growth phase were 
used, and the minimum acceptable R^ value for the regression line 
was 0.95. GT was then calculated by dividing 0.301 (equivalent to 
logio 2) by the slope of the line. LT was estimated from the time at 
which the regression line through the exponential part of the 
growth curve first exceeded the initial cell numbers after applying 
the hyperosmotic shift (i.e. the time taken to increase above 
starting numbers). We also determined the time taken for cells to 
resume exponential growth following the shift. This is referred to 
as 'adaptation time' (AT), and was estimated visually from the 
growth curves. To measure the amount of work to be done by 
E. coli to adjust to new conditions, relative lag time (RLT) was 
calculated by dividing LT or AT by GT [14,15]. 

Preparation of samples for transcriptomic and proteomic 
analyses 

In a separate set of hyperosmotic shift experiments, a rapid 
downshift of £. coli cultures from a^ 0.993 to a^ 0.967 was carried 
out as described above to prepare samples for cDNA microarray 
and 2D-LC/MS/MS analysis. An appropriate volume of cell 
cultures was harvested before osmotic upshift (i.e., control culture 
or non osmoticaUy-treated cells), and at 0 (i.e., immediately after 
the shift), 30, 80, and 310 min after the shift. Independent samples 
were then subjected to extraction of RNA, and both soluble and 
membrane proteins, according to the methods of Kocharunchitt 
et al. [9] . It also should be noted that the current proteomic study 
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used the same control samples (i.e., at 35°C a„ 0.993) as in tlie 
companion study [16]. The number of biological replicates 
performed for each time point in the transcriptomic and 
proteomic analyses is described in Table 1. 

Both transcriptomic (cDNA microarray) and proteomic (2D- 
LC/MS/MS) analyses were carried out as described by Kochar- 
unchitt et al. [9]. 

Microarray data analysis 

Gene expression was analyzed using GeneSpring GX 10.0 
software (Agilent Technologies, Palo Alto, CA). Array data was 
normalized using the Robust Multi-Array Analysis (RMA) 
algorithm [17]. Probesets with an intensity value in the lowest 
20th percentile among all the intensity values were excluded. The 
remaining probesets, with intensity values between the 20th and 
100th percentile, were used for downstream statistical analysis. 
Significance analysis was conducted using one way analysis of 
variance (ANOVA) using the multiple correction testing method of 
Benjamini and Hochberg [18] with a /"-value cut-off of <0.01. 
Genes were considered to be differentially expressed if the fold 
change was >2. 

MS/MS data analysis 

MS/MS data obtained from each protein sample of membrane 

and soluble fractions were processed by the Computational 

Proteomics Analysis System (CPAS), a web-based system built 

on the LabKey Server (v9.1, released 02.04.2009) [19]. The 

_ experimental mass spectra produced were subjected to a semi- 

I tryptic search against the combined databases of E. coli 0157:H7 

Sakai (5,318 entries in total) downloaded from the National Center 

for Biotechnology Information (NCBI, https:/ /www.ncbi.nlm.nih. 

gov/, downloaded 25.11.2008) using XITandem v2007.07.01 

[20]. These databases included the E. cofi 0157:H7 Sakai database 

(5230 entries, NC_002695.fasta) and two E. coli 0157:H7 Sakai 

plasmid databases, plasmid p0157 (85 entries, NC_002128.fasta) 

and plasmid pOSAKl (three entries, NC_002127.fasta). The 

£ o I parameters for the database search were as follows: mass tolerance 

- ^ for precursor and fragment ions: 10 ppm and 0.5 Da, respectively; 

fixed modification: cysteine cabamidomethylation (+57 Da); and 

._ ^ 8 "u ^ no variable modifications. The search results were then analyzed 

^ ^ al "° £ using the PeptideProphet and ProteinProphet algorithms from the 

i 2 1^ ^ -a Trans Proteomic Pipeline v3.4.2 [21,22]. All peptide and protein 

E S" c -i- c identifications were accepted at PeptideProphet and ProteinPro- 

^ '§ ^ c phet of &0.9, corresponding to a theoretical error rate of ^2% 

[21]. To assess a false-positive discovery rate for each dataset, the 

MS/MS spectra were searched against the target-reversed (decoy) 

database using the same search criteria as described above. The 

S _^ false-positive discovery rate of peptide identifications was then 

-§ .g estimated by dividing the number of spectra matching decoy 

1 p S peptides with the total number of spectra, and appeared to be less 

than 5% for all datasets (Table SI in File SI). 

3 c :! 5 All protein identifications that passed the above criteria were 

further assessed based on the confidence level of protein 

identifications across biological replicates of each fraction of E. 

c S S u <u S coli at each time point. The confidence level was based on the 

I 3 !5 -s ■£ ° 
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o number of unique peptides identified from one sample and the 
Q- S S 1 3 c number of replicates in which the protein was detected. 
S c c c ■- ^3- Specifically, proteins identified by more than one unique peptide 
= §2EgSc in at least one of the replicates were considered to have a 'high' 
^ ^ ^ ^ g ^ :2, confidence score. The 'intermediate' confidence level was assigned 
'° proteins with a single peptide hit that was detected in more than 
E "i £ "I — ^ ^ one rephcate. 'Low-confidence' proteins were considered to be 
z ^ < z % those identified by a single unique peptide and found in only 
one replicate. Only protein identifications with a 'high' and 
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Figure 1. Growth of £ co// 01 57:1-17 upon water activity downsliift. Growth response of £ coli was determined by viable count when 
subjected to a rapid downshift from a„ 0.993 to a„ 0.980 (■), 0.975 (□), 0.970 (•), 0.967 (O) and 0.960 (♦) at a constant temperature of 35°C. The 
time at which hyperosmotic shift was applied is indicated by a dotted line. Data points represent means ± standard deviations of at least two 
independent replicates. 
doi:l 0.1 371/journal.pone.0090422.g001 



'intermediate' confidence level (referred to as having a 'high 
confidence score') were accepted for further analysis. 

Following data filtering, all protein identifications with high 
confidence from membrane and soluble fractions of the same 
sample were combined to represent a 'total' proteome of the 
corresponding sample. This combined approach was carried out 
to improve reliability of the proteomic data. The values, as 
determined by linear regression analysis indicated a relatively 
stronger linear correlation between spectral counts (SpCs, the 
number of MS/MS spectra) for all possible pair-wise comparisons 
of replicates of pooled fractions (average = ().89±0.06) when 
compared to that of replicates of each fraction (average 
R^ = 0.64±0.20, and 0.85±0.06 for membrane and soluble 
fractions, respectively). Each of the total fractions across the 
biological replicates was then used to generate the final list of 
proteins considered to be representative of each time point of 
sampling. 

Protein abundance ratio and its significance 

Spectral count generated by 2D-LC/MS/MS analysis was used 
as a semi-quantitative measure of protein abundance [21]. The 
total abundance of each protein in a sample was obtained by 
pooling spectral counts detected for a given protein from 
membrane and soluble fractions of the same sample (i.e. the same 
biological replicate). Normalized spectral abundance factor 
(NSAF) was then applied to the pooled spectral count of each 
protein, according to ZybaUov et al. [23]. Briefly, the NSAF for a 
given protein is the pooled spectral count of the protein divided by 
its amino acid length (L), and normalized against the average 
SpC/L across all proteins in the corresponding dataset. The NSAF 
values were then averaged across all biological replicates at each 
time point of sampling to obtain the representative values of 
NSAF. 

Fold changes in protein abundance due to hyperosmotic shock 
were calculated as the logj ratio of the average NSAF value in the 
osmotically-treated sample over the average NSAF value of the 
control. To avoid errors during fold change calculation in cases 
where a zero value was obtained, a correction factor of 0.01 was 



added to each of the average NSAF values. Furthermore, 
difierences in protein abundance were statistically analyzed using 
the beta-binomial test implemented in R [24] . All proteins with a 
P- value of < 0.0 1 and at least 2 fold change were considered to be 

differentially abundant. 

Transcriptomic and proteomic data mining 

Information on identified genes and proteins such as protein 
and gene names, ECs numbers (locus tag), GI numbers, NCBI 
Reference Sequence (RefSeq), protein sizes and molecular masses 
were obtained from the UniProt knowledgebase (http:/ /www. 
uniprot.org/, accessed on 5.12.2010) and the NCBI database 
(accessed on 5.12.2010). Preliminary functions and properties of 
the genes and proteins were based on the EcoCyc database 
(http://www.ecocyc.org/, accessed on 13.12.2010), and the Kyoto 
Encyclopedia of Genes and Genomes (KEGG) database (http:// 
www.genome.jp/kegg/, accessed on accessed on 4.12.2010). 
Protein or gene names in conjunction with ECs numbers were 
used as a unique identifier for genes and/ or proteins. 

Predefined set enrichment analysis 

To determine changes in overall expression of a predefined set 
of genes or proteins due to osmotic upshift, a t-test based method 
of Boorsma et al. [25], known as T-profiler, was performed on the 
log2 expression ratios. The predefined sets of genes or proteins 
used for T-profiler analysis were based on the functional role 
categories and/ or metabolic pathways of the following databases: 
the JCVI Comprehensive Microbial Resource (JCVI CMR) 
database (http:// cmrjcvi.org/ cgi-bin/CMR/CmrHomePage.cgi, 
accessed on 5.12.2010), the EcoCyc database collection (http:// 
www.ecocyc.org/, accessed on 13.12.2010), and the KEGG 
database (http://www.genomejp/kegg/, accessed on 4.12.2010), 
as well as based on the lists of genes or proteins previously known 
to be positively controlled by major regulons, including CpxRA 
[26], RpoE [27], RpoH [28], RpoS [29-31], and the Res 
phosphorelay system [32]. The T-value obtained from the analysis 
was determined only for sets that contained at least five members, 
and its significance was established by using the associated 
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Table 2. Summary of the growth kinetics of exponential phase £ coli 0157:H7 Sakai subjected to a sudden downshift in water 
activity from a^ 0.993 to aw 0.980, 0.975, 0.970, 0.967 and 0.960, as determined by viable count. 





Relative lag time (RLT) 


Water activity 
downshift 


Generation time, 
GT (h) 


Lag time, LT (h)* 


Adaptation time, 
AT (h)" 


LT/GT 


AT/GT 


a„ 0.980 


0.45 


0.69 


0.68 


1.53 


1.51 


a„ 0.975 


0.67 


1.55 


0.94 


2.31 


1.40 


a„ 0.970 


0.73 


2.04 


1.27 


2.79 


1.74 


a„ 0.967 


1.41 


4.74 


1.53 


3.36 


1.09 


a„ 0.960 


1.68 


10.57 


2.13 


6.29 


1.27 



^Lag time is calculated by linear regression as the time tal<en to increase above starting numbers. 
"^Adaptation time is estimated visually as the time taken for cells to resume exponential growth. 
doi:l 0.1 371 /journal.pone.0090422.t002 



two-tailed P-value. All predefined sets with a P-value less than 0. 1 
were considered to be statistically significant. 

Clustering analysis of transcriptomic and proteomic data 

To determine similarities in transcriptomic and proteomic 
profiles of E. coli during exposure to a sudden downshift in a^, 
hierarchical clustering analysis (HCA) was performed on the T- 
values calculated for JCVI CMR functional categories. The T- 
profiler results for the transcriptome and proteome oiE. coli during 
steady-state growth at 25°C a„ 0.967 (data of Kocharunchitt et al. 
[9]) were also included in the clustering analysis for comparison. 
The degree of similarity was calculated with the Euclidean 
distance metric and complete linkage as a clustering method, using 
Cluster 3.0 software [33]. The clustering result was then visualized 
in Tree View version 1.1.3 [34] . 

Transcriptomic and proteomic data accession numbers 

The transcriptomic data are available in the ArrayExpress 
database (www.ebi.ac.uk/arrayexpress) [35] under accession 
number E-MTAB-2011. The mass spectrometry proteomic data 
are deposited to the ProteomeXchange Consortium (http:// 
proteomecentral.proteomexchange.org) via the PRIDE partner 
repository [36] with the dataset identifiers PXD000.583 and 
PXD000541. 

Results and Discussion 

Growth of E. coli 0157:H7 Sakai upon a„ downshift 

Growth of exponential phase E. coli (ODgoo of ~0.1±0.01) was 
assessed in a nutrient rich medium when subjected to a rapid 
downshift from a„ 0.993 to a„ 0.980, 0.975, 0.970, 0.967 and 
0.960 at a constant temperature of 35°C. All downshifts resulted in 
a growth interruption before cells resumed their growth at a slower 
rate (Figure 1). Estimates of kinetic parameters for each a„ 
downshift are summarized in Table 2. The generation times 
obtained were consistent with those expected from the predictive 
model of Ross et al. [37], and both generation and lag times 
increased with respect to the magnitude of the osmotic shift. The 
RLTlx/gt response of E. coli increased in a linear manner 
(y= -225.64x-l.84, R2 = 0.89) with the magnitude of the shift. 
This is consistent with previous studies, reporting that cells shifted 
to progressively harsher conditions require more adaptation work 
to be done, i.e. RLT become larger [13,38,39]. 

E. coli appeared to have a complex pattern of growth in response 
to sudden downshifts in a^. A rapid, but transient decrease in 
bacterial numbers was observed immediately after a shift to lower 



a„ (Figure 1). This indicates that cells might lose their culturability 
on enumeration media, possibly due to injury. The degree of loss 
of culturability also became progressively larger with the larger 
shifts in a„, although this loss of culturability occurred at a similar 
rate for all a„ tested. However, the loss of culturability was 
reversed within the lag period at a rate faster than that 
characteristic of exponential growth under the given conditions 
(Figure 1). Furthermore, it was evident that during treatment at a„ 
0.967 and 0.960, the exponential growth occurred before the 'lag 
period' (when defined as the time taken to increase above starting 
numbers) was resolved. This suggests that cell death occurred 
upon the shift at the two lowest a„ tested. A similar injury and 
recovery phenomenon has been reported in other studies following 
the application of sub-lethal stressful treatments [40-42] . Howev- 
er, our results are in contrast to the study of Mellefont et al. [38] in 
which abrupt osmotic shifts of Salmonella Typhimurium caused 
inactivation (rather than loss of culturability) of a sub-population 
followed by renewed exponential growth. 

The complexity of the response of £. coli to osmotic upshifts (i.e. 
irregularly shaped growth curves arising from the initial decline 
and subsequent recovery in cell numbers) may contribute to 
unusual estimates of RLT, although a consistent and reproducible 
pattern of the RLTlx/gt response was obtained in this study. 
Death of cells would mask the time when the lag period is actually 
resolved, resulting in overestimation of lag time. This led us to 
determine adaptation time (AT), i.e. the time taken for bacterial 
cells to resume exponential growth after hyperosmotic shift. The 
adaptation times increased with the magnitude of the osmotic shift, 
a similar trend to that observed in lag times (Table 2). However, 
when RLT estimates were calculated using the adaptation time, 
there was no systematic variation in RLTsat/gt over the a„ range 
tested (y= 14.53x+1.73, R^ = 0.20). This implies that surviving E. 
coli cells require a similar amount of work to be done to adapt to a 
given decrease in a„ and, thus, is inconsistent with the observed 
pattern of RLTlx/gt response and previous studies in which shifts 
to progressively harsher environments would require more 
adaptation work to be done [13,38,39]. Further investigations at 
more extreme osmotic upshifts together with validation of our data 
by culture-independent methods are necessary to confirm the 
present observations. 

Molecular response of £ coli 0157:H7 Sakai to 
hyperosmotic shock 

To understand better the potential mechanisms adopted by 
E. coli to respond to osmotic shock and hyperosmotic stress, the 
present study employed transcriptomic and proteomic analyses in 
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parallel in independent experiments to examine temporal alter- 
ations in mRNA expression and protein production in E. coli upon 
a rapid downshift from a^^ 0.993 to a^ 0.967 at a constant 
temperature of 35°C. Bacterial cells were harvested before 
hyperosmotic shift, and 0 (i.e. immediately after the shift), 30, 
80, and 310 min after the shift (Figure 1). It also should be noted 
that these time points were chosen with the aim to characterize the 
physiology of E. coli during dynamic changes in growth kinetics 
induced by hyperosmotic shift. Specifically, the samples taken at 
time 0 and 30 min were obtained during the period in which the 
loss of cell culturability was observed, whereas the samples at time 
80 and 3 1 0 min respectively reflected the physiological state of E. 
coli during the 'recovery' period and during growth after the shift. 
Independent samples were processed for mRNA, and membrane 
and soluble protein extractions as appropriate for cDNA 
microarray and 2D-LC-MS/MS analyses, respectively. By using 
multidimensional LC /MS/MS analysis together with data filter- 
ing, protein identifications with high confidence were obtained for 
total fractions of £. coli at each time point (Table 1). 

Changes in gene expression and protein abundance due to 
hyperosmotic stress were determined by comparing the profiles 
obtained after the shift to that of a reference sample taken prior to 
the shift. Table 1 describes the number of genes and proteins with 
significant changes in expression level at each time point. It was 
found that the number of up-regulated elements (referred to here 
as genes and proteins) was dominant o\ er those with decreased 
expression level. A complete list of these- genes and proteins with 
their comprehensive annotations, as well as their comparison to 
previously pubhshed studies is shown in Table S2 in File SI. 

Exposure of E. coli to a rapid downshift in a,^ evokes a highly 
complex regulatory process involving expression changes of 
different groups of genes and proteins, which were condensed 
using the T-profiler analysis (see Tables SS-S.^ in File SI). These 
genes and proteins are involved in a number of functional 
categories and metabolic pathways, as well as play an important 
role in the stress responses that are controlled by regulons, as 
discussed below. 

Comparison of transcriptome and proteome 

To determine the level of correlation between the transcrip- 
tomic and proteomic data, the Ust of differentially expressed genes 
and proteins at each time point was compared. It was found that 
the level of correlation was low with an average of only 
4.32% ±7.09% of genes and proteins matched across all time 
points (Table 1). This was despite that a number of trends were still 
observed in the data (discussed below). Consistent with these 
results, several studies have reported poor correlation between 
gene and protein expression profiles [9,43,44]. Several possible 
explanations for this have also been suggested, including different 
controlling mechanisms at the transcriptome and proteome levels, 
post-transcriptional mechanisms affecting the protein translation 
rate, differences occurring between the half-lives of specific 
proteins or mRNAs, and the intracellular location and molecular 
association of the protein products of expressed genes [45] . 

To further determine the similarity between gene and protein 
profiles across different time points, the clustering analysis was 
performed on the T-profiler results obtained in the present study 
(i.e. both transcriptomic and proteomic data) and that of 
Kocharunchitt et al. [9] on the transcriptomic and proteomic 
responses of E. coli during exponential growth at a similar osmotic 
stress condition (25°C a^ 0.967) (Figure 2). It was found that gene 
and protein profiles could be grouped into three major clusters. 
The first cluster is comprised of the gene profile at 0 min after 
downshift in SLy, and the protein profiles within the first 80 min 



following the shift, whereas the second cluster contains the gene 
profiles at time 30 and 80 min, and the protein profile at time 
310 min and that of steady-state growth at 25°C a^ 0.967. For the 
last cluster, it consists of the gene profile at 310 min and that of 
steady-state growth at 2.^°C a„ 0.967. These findings, taken 
together, indicate that the transcriptomic data at a given time 
point was a reflection of subsequent, rather than concomitant 
proteomic data, providing a valid explanation for the low level of 
correlation between gene and protein profiles observed here (see 
above). The apparent alignment of a given transcriptomic data 
with proteomic data at a later time point has also led to the 
suggestion that E. coli cells might take longer to adjust their gene 
expression (i.e. transcription) in response to hyperosmotic shock 
and that change in the processes of translation and post- 
translational modification might be the potential mechanism used 
by bacteria to quickly adapt or respond to such environmental 
change. In keeping with this, Rossouw et al. [46] have reported a 
similar observation for the transcriptome and proteome of 
Saccharomyces cerevisiae during alcohoUc fermentation. Furthermore, 
it should be noted that the last two clusters imply that the 
transcriptome and proteome of E. coli at 310 min after the shift 
were similar to those of E. coli grown at high osmolality (25°C a^^, 
0.967). This supports earlier observations that E. coli had already 
estabhshed balanced growth typical of growth under such stress. 

Despite the discordant alignment of transcript and protein levels 
observed here, previous reports have indicated the need for 
complementary profiling of transcriptome and proteome to fully 
develop an understanding of the functional architecture of 
genomes and gene networks [47,48]. 

Transcription and translation 

Bacteria generally control the cellular content of transcriptional 
and translational machineries according to their physiological 
state. Faster-growing cells synthesize proteins more quickly, and 
the levels of ribosomes and ribosome-associated proteins correlates 
well with growth rate [49,50] . Reduced growth rate as a result of 
environmental conditions including high osmoIarUty has an effect 
on gene expression and protein production [9,51]. However, the 
present findings did not show a strong alteration in expression of 
genes and proteins associated with RNA polymerase in response to 
Zy, downshift (Table S4 in File SI). This was despite the 
transcriptomic analysis revealing a significant down-regulation of 
pcnB encoding poly(A) polymerase I (1.06 logg-fold) only at 80 min 
of the treatment. On the other hand, it was found that the level of 
translational machinery was affected after prolonged exposure to 
osmotic stress. Several ribosomal proteins exhibited a decrease in 
their abundance at 3 1 0 min after the treatment, as indicated by a 
significant negative T-value of 4.39 (Table S4 in File SI). This 
response agrees well with the observations of Kocharunchitt et al. 
[9] who reported a reduced translational capacity in osmotically- 
adapted cells. 

Carbohydrate catabolism and energy generation 

A sudden downshift in a,^ was genercdly found to have no major 
effect on expression of genes and proteins involved in carbohy- 
drate metabolism during the adaptation period following hyper- 
osmotic shock (i.e. at time 0, 30 and 80 min after the shift). Only 
three genes encoding the enzymes in the glycolysis/glucxmeogen- 
esis pathway exhibited a significant change in transcription level, 
including JbaB (4.60 logj-fold increase), tpiA (1.03 log2-fold 
decrease) and jhhA (1.15 log2-fold decrease). On the other hand, 
when the osmotically-treated cells had resumed exponential 
growth (i.e. at time 310 min), several genes and proteins of 
the glycolysis and gluconeogenesis pathways appeared to be 
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Figure 2. Hierarchical clustering analysis of transcriptomic and proteomic data. The heat map represents the l values calculated for 
transcriptomic (T) and proteomic (P) profiles based on JCVI CMR functional categories (row), and is linked by a dendrogram representing clustering of 
these profiles obtained from £ coli upon hyperosmotic shift, and those during exponential growth at 25°C a„ 0.967 [9] (top). The color code is as 
follows: red indicates positive 7-values; green indicates negative T-values; and grey indicates not determined. 
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up-regulated, although significant changes in their overall 
expression were not evident in the T-profiler analysis (Table S4 
in File SI). These included JbaB/FhaAB, pfkB/FfkAB, mdh/Mdh 
and maei?/MaeB for which the expression level increased by up to 
6.40 log2-fold. Similarly, up-regulation of a number of genes and 
proteins involved in the pentose phosphate pathway was evident at 



this time point. Among these, transaldolase A (talA/TalA) and 
transketolase B [tktB/TktS), which are involved in the non- 
oxidative branch of the pentose phosphate pathway, exhibited a 
significant up-regulation. Regulation of these genes is also known 
to be positively controlled by RpoS at the transcriptional level 
(see below) [31]. 
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E. coli typically utilizes three alternative modes for energy 
generation, namely aerobic respiration, anaerobic respiration and 
fermentation [52]. It has been shown that the aerobic respiration 
of E. coli upon exposure to hyperosmotic shock is severely inhibited 
as a consequence of plasmolysis, and that the respiratory activity is 
recovered during d(;plasmolysis [53]. In the present study, the 
transcriptomic and proteomic data were inconsistent with regard 
to responses related to oxidative phosphorylation (Table S4 in File 
SI). A number of genes encoding the enzymatic complexes 
involved in oxidative phosphorylation were down-regulated at 
80 min of the tr(;atmc'nt, as indicated 1)y a significant negative T- 
value of 2.43. These included genes nuoKLMM (1.23 log2-fold 
decrease on average) encoding the membrane components of 
NADHiubiquinone oxidoreductase I. By contrast, the proteomic 
analysis indicated a significant increase in overall abundance of 
proteins involved in oxidative phosphorylation within the first 
80 min. A common set of proteins observed at these time points 
were the components of NADH:ubiquinone oxidoreductase I 
(NuoCFI; 2.94 logg-fold on average), succinate dehydrogenase 
(SdhA; 2.77 loga-fold on average), and ATP synthase (AtpB; 3.32 
log2-fold on average). These inconsistent observations from two 
different methodologies may be in line with those of Gygi et al. 
[45] who compared mRNA and corresponding protein levels, as 
discussed earlier. However, E. coli appeared to increase aerobic 
respiration after prolonged exposure to hyperosmotic stress, as 
indicated by significant positive T^values for the TCA cycle and 
glyoxylate bypass pathways (Table S4 in File SI). The transcrip- 
tomic study revealed a significant increase in expression of several 
genes [oceAB, acriA, sdliACD and gltA; 3.49 log2-fold on average) 
encoding the enzymes in the TCA cycle and glyoxylate bypass at 
80 and 310 min of the stress, whereas almost all the TCA cycle 
and glyoxylate bypass-related proteins were found to be signifi- 
candy up-regulated only at time 310 min in the proteomic analysis 
(Table S2 in File SI). 

Taken together, the increased expression of genes and/or 
proteins involved in the major processes of carbohydrate 
catabolism and energy generation at 310 min after the downshift 
may be to compensate for the reduction in activity of these 
processes at low a,^, or may reflect a high level of energy 
production together with an increase in level of precursors for the 
biosynthesis of various macromolecules. The latter corresponds 
well with the observation that E. coli resumed their growth at or 
before this time point (Figure 1). Our previous study also indicated 
an increase in metabolic activities to meet an energy requirement 
during exponential growth at high osmolality (25°C a„ 0.967) [9]. 

Two-component regulatory system 

The mechanisms involved in osmosensing and subsequent 
regulatory responses in E. coli are strongly controlled by two- 
component regulatory systems [6,54]. The EnvZ/OmpR system 
regulates the transcription of genes encoding porins OmpC and 
OmpF in response to changes in medium osmolality. This 
regulation is characterized by an altered ratio of these porins, 
whereby the synthesis of OmpC increases, and that OmpF 
decreases under hyperosmotic stress [55,56]. However, we did not 
observe a strong activation of this system (Table S2 in File SI), 
although a significant increase in the abundance of OmpC and 
OmpR was observed at time 0 and 310 min of the osmotic stress, 
respectively. 

Amino acids biosynthesis 

Biosynthetic genes and proteins for various amino acids had 
altered expression levels after an abrupt downshift in a^^ (Tables S2 
and S3 in File SI). Among these, the mRNA level of AisCDF, 



which encode key enzymes involved in histidine biosynthesis, was 
significandy reduced at 80 min of the shift (up to 1.14 logj-fold). 
This was accompanied by the observation based on the T-profiler 
analysis that a significant decrease in the overall expression of 
histidine biosynthetic genes occurred at this time (Table S3 in File 
SI). The biosynthesis of histidine is known to be an energetically 
expensive process, requiring 41 ATP molecules per histidine 
molecule made [57]. Previous studies have also demonstrated that 
constitutive expression of the his operon caused a number of 
physiological changes in E. coli, including growth inhibition in 
high-salt media [58,59]. Therefore, the apparent decrease in the 
activity of the histidine biosynthetic pathway might imply that E. 
coli suppresses this biosynthetic pathway to prevent an inappro- 
priate use of the cellular energy, possibly for re-aUocating the 
energy to repair cell damage as observed during the recovery 
period (i.e. at time 80 min) (Figure 1). 

On the other hand, the present study revealed an induced 
activity of the methionine biosynthetic pathway in osmoticalh - 
adapted cells (i.e. at 310 min of the treatment), as indicated by a 
significant up-regulation of melEH genes and MetBEH proteins (an 
increase in abundance of methionine transporters was also 
observed, see below). Methionine is required for initiation and 
elongation of proteins, biosynthesis of polyamines, purines and 
pyrimidines, and various methylation reactions [60,61]. A 
previous study has also demonstrated that the growth of E. coli 
at elevated temperatures is limited by the availability of 
endogenous methionine [60]. This has led to the suggestion that 
an increase in the cellular level of methionine might contribute to 
the maintenance of balanced growth rate of E. coli at high 
osmolality. However, Kocharunchitt et al. [9] did not find up- 
regulation of the methionine biosynthetic pathway during growth 
at a similar osmotic stress condition of 25°C Hy, 0.967. This might 
be due to the diflFerence in the conditions tested in this study (i.e. 
35°C 0.967). Furthermore, the T-profiler analysis r('\ealed a 
significant increase in overall expression of genes and proteins 
responsible for the biosynthesis of amino acids in the pyruvate 
family at the time point at which osmoticaUy-treated cells resumed 
their growth after osmotic upshift (i.e. at time 310 min) (Table S3 
in File SI). Among these, branched-chain amino acids biosynthetic 
pathway genes and proteins (zfoiiDG/IlvBDG) exhibited a signif- 
icant up-regulation. This agrees well with the observation on the 
adaptive response of E. coli to osmotic stress (25°C 0.967) [9]. 
Taken together, the apparent increase in expression of several 
genes and proteins related to the amino acid biosynthetic pathways 
(i.e. methionine, isoleucine and valine) supports the previous 
hypothesis that an increase in intracellular contents of these amino 
acids may be part of the responses to enhance growth of E. coli 
under osmotic stress [9]. Horinouchi et al. [62] have also shown 
that supplementing culture media with specific amino acids (i.e. 
tryptophan, histidine and branched-chain amino acids) promoted 
the growth of E. coli under another type of stress (i.e. ethanol 
stress). 

Transport functions 

Accumulation of compatible solutes through uptake from the 
external environment by active transport is one of the most 
important responses for bacteria to counter hyperosmotic stress 
[7]. Indeed, the present study found a significant increase in 
activity of a number of transporters for compatible solutes upon 
exposure to osmotic upshift (Table S2 in File SI). Although the 
high-affinity potassium ion (K"*) transport system (Kdp) did not 
change in its level after the shift, abundance of the low-affinity 
system Kup, was up-regulated by up to 3.57 log2-fold at aU time 
points. The potassium channel protein Kch also exhibited an 
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increase in its abundance at time 80 min and 310 min (2.55 log2 
-fold on average). Correspondingly, the expression of genes (gltjl) 
encoding transporters for glutamate (the potassium ion counter 
ion) was found to be induced at these time points (1.47 log2-fold on 
average). Other transport systems that were up-regulated in both 
transcriptomic and proteomic analyses included those responsible 
for the uptake of glycine betaine/proline (proVWX operon; at 30, 
80 and 310 min of the treatment) and other osmoprotectants 
(yehWXrZ operon; at 80 and 310 min). In addition, NhaAB 
sodium ion/proton antiporters exhibited a significant increase in 
their abundance throughout the course of osmotic stress. These 
antiporters are known to be impKcated not only in osmotic 
adaptation, but also pH homeostasis [63]. 

Apart from the expression of genes and proteins for the uptake 
of several compatible solutes as already mentioned, exposure of i?. 
coli to a sudden downshift in a„ also caused a significant alteration 
in the activity of other transport systems. It was observed that 
abundance of several transport proteins for phosphate and amino 
acid was induced. These included the proteins involved in the 
uptake of arginine (ArtQMP; at 310 min of the treatment) and 
methionine (MetQN; at all time points). For the group 'peptide 
and nickel transporters', production of SapCF, which are 
components of the SapABCDF uptake system for peptides, was 
induced at 30 and 310 min after the a„ downshift. Proteins 
involved in the transport of murein oligopeptides (OppBD) also 
exhibited a significant increase in expression at time 310 min. 

The transcriptomic analysis revealed a significant down- 
regulation of several genes, encoding transporters for minerals, 
and inorganic and organic ions, such as sulfate/thiosulfate (cysWP, 
at 80 min after a„ downshift) and molybdate [modABC; at 
310 min). Transcription of the potAB operon involved in the 
uptake of polyamines, putrescine and spermidine was also reduced 
at 80 and 310 min following the shift. This may indicate a 
decrease in the cellular level of polyamines. Lack of polyamines 
has previously bc(;n slio\\n to be linked to abnormal growth and 
oxidative stress-induced damage [64]. Furthermore, genes encod- 
ing monosaccharide and oligosaccharide transporters were found 
to be up-regulated. These included those responsible for the 
uptake of xylose (xylFHG; at 80 min) and maltose [malEGFK; at 
310 min). It has previously been reported that catabolism of 
maltose can contribute indirecdy to enhance the levels of 
endogenously synthesized osmolytes such as glutamate and jV- 
acetylglutaminylglutamine amide, facilitating growth at elevated 
osmolarities [65]. 

Increased production of several proteins associated with the 
uptake system for metallic cations, iron-siderophores and vitamin 
B12 was observed in response to osmotic upshift. TonB and 
ExbBD, which make up the TonB-ExbBD energy transducing 
system to provide the energy source required for the uptake of 
iron-siderophore complexes and vitamin B12 across the outer 
membrane, were significantly up-regulated throughout the treat- 
ment. This was accompanied by the obscr\ atioiis that FhuCD 
were significantiy up-regulated (4.42 log2-fold on average) at all 
time points, indicating increased activity of the TonB-dependent 
iron (III) hydroxamate transport system. Abundance of FepA, a 
subunit of the TonB-dependent uptake system for ferric enter- 
obactin also appeared to be significantly induced (2.63 log2-fold on 
average) within the first 30 min of the shift. In keeping with these 
findings, the response of E. coli during steady-state growth at a 
similar stress condition (25°C a„ 0.967) revealed an increase in 
overall abundance of proteins involved in the transport of cations 
and iron compounds [9]. The enhancement of iron uptake 
observed here may be associated with an increase in levels of 
intracellular hydrogen peroxide (H2O2) and superoxide (O2"), 



which invoke the oxidative stress response [62]. This also 
corresponds well with the present observations in which a number 
of genes and proteins involved in defense mechanisms against 

oxidative stress were significantly up-regulated by osmotic stress 
(see below). In addition, BtuE a component of the TonB- 
dependent vitamin B12 transport system exhibited a significant 
increase in its abundance of up to 2.45 log2-fold within the first 
80 min of the stress. 

Fatty acids and lipids metabolism 

E. coli membranes typically comprise three major phospholipids, 
including phosphatidylethanolamine (PE), phosphatidylglycerol 
(PG) and cardiolipin (CL). Growth at low d^^ typically leads to 
changes in membrane lipid composition by increasing the 
proportion of anionic (PG and CL) over zwitterionic (PE) 
phospholipids. Such changes are to prevent the membrane lipids 
from adopting a non-bUayer phase, which otherwise could disrupt 
membrane function [66]. However, neither transcriptomic nor 
proteomic data provided strong evidence supporting this adaptive 
strategy. A significant increase in overall abundance of enzymes in 
the biosynthetic pathway of phospholipids was observed within the 
first 30 min of imposition of the osmotic stress (Table S4 in File 
SI). Phosphatidylserine synthase (PssA) and cardiolipin synthase 
(Cls) were found to be significantiy up-regulated (2.94 log2-fold on 
average). PssA and Cls proteins are responsible for the de novo 
synthesis of PE and CL, respectively, suggesting that increased 
production of these phospholipids occurs. Despite that the 
physiological role of PssA under osmotic stress is unclear, PE has 
previously been shown to be important for acid adaptation [67]. 

Hyperosmotic stress was found to have no significant effects on 
overall expression of genes and proteins involved in the major 
pathways for fatty acid biosynthesis (Table S4 in File SI). This was 
despite the observation that the abundance of FabZ was 
significantiy induced by up to 4.03 log2-fold at all time points of 
the stress. The FabZ enzyme plays an important role in the 
biosynthetic pathways of both saturated and unsaturated fatty 
acids. In contrast, the transcriptomic analysis revealed a significant 
up-regulation of the cfa gene, which encodes cyclopropane fatty 
acyl phospholipid synthase at time 80 and 310 min. The increase 
in cyclopropane fatty acid content in the cell membrane has 
previously been demonstrated to assist cells to maintain intracel- 
lular pH homeostasis by reducing membrane permeability to 
protons. This lipid modification provides protection against acid 
stress and other stress conditions such as high salt concentration 
and ethanol [68-70]. Furthermore, our observations on fatty acid 
composition are consistent with the study of GuiUot et al. [7 1] on 
the response of Lactococcus lactis to osmotic stress, in which that 
organism was reported to increase the level of cyclopropane fatty 
acids, whereas the unsaturated-to-saturated fatty acids ratio 
remains unchanged. 

Several genes, encoding the key enzymes for lipopolysaccharide 
biosynthesis (i.e. IpxH and waaA) were down-regulated with a 
significant negative T-value at 80 and 310 min after hyperosmotic 
shift (Table 84 in File SI). The apparent reduction of lipopoly- 
saccharide biosynthesis might indicate that outer membrane 
instability occurs during adaptation to hyperosmotic stress. 
Consistent with this, previous studies have reported that a defect 
in lipopolysaccharide biosynthesis leads to the lack of a continuous 
lipopolysaccharide layer in the outer membrane, causing increased 
susceptibility of bacterial cells to hydrophobic antibiotics [72,73]. 
Furthermore, it has been demonstrated that the envelope stress 
caused by the defective biosynthesis of lipopolysaccharide increas- 
es the biosynthesis of the exopolysaccharide colanic acid [74,75]. 
This, indeed, agrees well with earher observations, indicating that 
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the Res system-regulated colanic acid biosynthesis becomes 
activated (see below). 

Cell structure components 

In response to hyperosmotic stress, E. coli increased expression 
of several genes and proteins involved in the Res phosphorelay 
system that regulates the biosynthesis of colanic acid [32] (Table 
S5 in File SI). The T-profiler results revealed that a significant 
increase in overall expression of genes known to be induced by Res 
regulon [32] occurred from 30 min of the osmotic treatment 
onward, whereas several Res-dependent proteins were significandy 
up-regulated only at the time point at which E. coli had resumed 
growth (i.e. at time 310 min). These genes and proteins were also 
found to be amongst the most highly up-regulated in the present 
study (Table S2 in File SI). Consistent with these findings, 
Kocharunchitt et al. [9] also demonstrated strong up-regulation of 
Res-dependent elements together with a high level of colanic acid 
production in E. coli cells during steady-state growth under a 
similar stress condition (25°C 0.967). The importance of colanic 
acid has frequendy been described as protecting cells against a 
variety of stresses, including osmotic stress [76-78], and has been 
shown to be involved in biofilm formation [79]. Although the 
physiological role of colanic acid is not well understood, it is 
thought that colanic acid expressed on cell surfaces simply 
provides a physical barrier to protect ceUs from hostile environ- 
ments [78,80]. Allen et al. [81] has reported that colanic acid 
confers a strong negative charge to the cell surface. This negatively 
charged cell surface has led to the suggestion that colanic acid may 
help E. coli to maintain hydration of the cell surface, and to 
preserve the membrane hpids in the proper bilayer phase, as part 
of the adaptive strategies in response to such stress [65,66]. 
However, the findings of Kocharunchitt et al. [9] indicated that 
colanic acid is not required for growth and survival under osmotic 
stress. A comprehensive understanding on the function of colanic: 
acid is, therefore, needed to understand the indutxd activity of the 
Res system-controlled colanic acid biosynthesis under osmotic 
stress. 

Virulence properties 

Adaptation of E. coli to high osmolality appeared to be 
associated with increased virulence. It was evident in the 
transcriptomic analysis that several genes encoding virulence 
factors in the locus of enterocyte eflfacement (LEE) pathogenicity 
island (PAI) were significantly up-regulated. These included kr 
encoding a regulator required for activation of LEE (from 80 min 
of the stress onward), eae encoding the outer membrane adhesin 
intimin (only at time 80 min), tir encoding the translocated intimin 
receptor (from 30 min onward), and esc operon encoding major 
components of type III secretion system (only at time 310 min). In 
contrast, the proteomic data did not provide strong evidence of 
increased abundance of virulence proteins (Table S2 in File SI). 
The obser\'ed increase in transcript levels of virulence genes, 
however, supports the idea that bacteria typically enhance their 
virulence after exposure to environmental stresses [6]. 

Bacterial chemotaxis and motility 

An immediate increase in overall expression of gc-ncs and 
proteins involved in chemotaxis and motility was observed after 
osmotic upshift, as indicated by significant positive T-values (Table 
S3 in File SI). This was followed by a general decrease in their 
level over the course of osmotic treatment. Accordingly, our 
transcriptomic data revealed a significant down-regulation of rpoF, 
which encodes a sigma factor responsible for controlling the 
transcription of genes involved in chemotaxis and motility [6], at 



80 and 310 min of the treatment. It was also found that 
chemotaxis and motility genes and proteins were amongst the 
most heavily down-regulated in osmotically-adapted cells (i.e. at 
time 310 min) (Table S2 in File SI). Kim et al. [82] has also 
reported a transient increase in expression of flagella biosynthetic 
elements in response to osmotic shock. This induced activity may 
be an initial response for cells to move away from osmotic stress 
[83] . Furthermore, the strong down-regulation of chemotaxis and 
motility systems observed here supports recent suggestions that 
these systems are the most dispensable functions during steady- 
state growth at high osmolality (25°C a„- 0.967) [9]. This adaptive 
response of E. coli is in keeping with the fact that expression of 
genes encoding flagella are negatively controUed by the master 
regulator RpoS [29], and the Res system-activated colanic acid 
biosynthesis [84]. 

General stress response 

Exposure of E. coli to a sudden downshift in a„ resulted in 
induced activity of the master stress regulator RpoS. This was 
evident in the transcriptomic data in which the mRNA level oirpoS 
significandy increased up to 3.40 logv-fold from 30 min of the shift 
onward, whereas the proteomic analysis revealed a significant up- 
regulation of RpoS only at time 80 and 310 min (4.66 log2-fold on 
average). Accordingly, the T-profiler analysis of both datasets 
showed a large increase in expression of several genes and proteins 
previously known to be induced by the RpoS regulon [29-3 1] over 
the course of treatment (Table S5 in File SI). Many of these have 
been described to be involved in osmoprotection [4,51], including 
osmotically inducible proteins {asm operon; at 30, 80 and 310 min 
of the treatment); (^j/Dps, which is responsible for the protection 
of DNA (at 80 and 310 min); otoAS/ OtsAB, which are involved in 
the de novo synthesis of trehalose (at 80 and 310 min); and a group 
of elements with undefined functions (f/a_B/ElaB, jaiiiVYciEF, 
j^fliW YgaM and j^TiJ/YjbJ; at 80 and 310 min). The transcrip- 
tomic study also found that a number of RpoS-dependent genes, 
which are known to be associated with defense mechanisms 
against other stresses [31,85], increased expression at 80 and 
310 min of the treatment. These included genes within the acid 
resistance fitness island (gadABWX and hdeBD), and those involved 
in response to oxidative damage [katE and mJABD). The proteomic 
analysis, however, revealed up-regulation of proteins involved in 
acid (GadAB; 2.47 log2-fold on average) and oxidative stress 
responses (KatE and SufAD; 5.35 log2-fold on average) only at 
time 310 min. Taken together, these observations support 
previous reports, demonstrating that the general response network 
estabhshed by the RpoS regulon typically provides cross-protec- 
tion against diverse stress conditions [9,29-31]. 

Apart from activation of the RpoS regulon as described above, 
E. coli cells transientiy induced additional stress responses or those 
elements controlled by other regulons (i.e. the RpoE and RpoH 
regulatory systems) during adaptation to hyperosmotic stress 
(Table S2 in File SI). A significant up-regulation of rpoE/'SjpoE 
was observed at 30 min (3.19 log2-fold; in the proteomic analysis) 
and 80 min (2.45 log2-fold; in the transcriptomic analysis) of the 
osmotic stress, while ?^o// expression significandy increased (1.81 
log2-fold on average) at time 30 and 80 min. The RpoE and 
RpoH regulons are known to be responsible for the control of 
protein misfolding in the cell envelope and cytoplasm, respectively, 
and it has been reported that the transcription of rpoH is RpoE- 
dependent [28,86]. Furthermore, the majority of genes whose 
expression is under the control of RpoH were found to be 
significandy up-regulated during the recovery period after 
hyperosmotic shock (i.e. at time 80 min). These included those 
encoding molecular chaperones (dnaK, htpG and ibpAB; 2.51 log2 
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-fold on average), which function in facilitating protein re-folding 
and stabilization, and preventing aggregation of misfolded 
proteins. DnaK is a member of the major DnaK/DnaJ/GrpE 
chaperone system, and has previously been reported to be 
involved in adaptation to osmotic and oxidative stresses [52,87]. 
Although other chaperone proteins (i.e. HtpG and IbpAB) have 
only minor activity in E. coli [88], previous studies have suggested 
that they play a significant role in de novo protein folding under 
stressful conditions, presumably by expanding the activity of the 
major system to interact with newly synthesized polypeptides 
[88,89]. 

The apparent response of £. coli cells to osmotic upshift through 
a transient induction of the RpoE and RpoH systems, and 
activation of the RpoS regulon reflects their robust ability to 
activate multiple defense mechanisms, and may also underlie the 
observation of a complex pattern of microbial growth behavior in 
response to such stress. In support of this, it has previously been 
suggested that increased activity of the RpoE and RpoH regulons 
may serve as a multipurpose emergency response to repair protein 
misfolding and facilitate the proper folding of newly synthesized 
proteins after loss of turgor (as would be expected during the injury 
and recovery period), while simultaneously activating the master 
stress regulator RpoS to mediate long-term adaptation at high 
osmolality [9,89]. 

The present study pro\ides evidence suggesting that E. coli 
might form filaments during adaptation to an abrupt downshift in 
a„. It was observed in the transcriptomic data that expression of 
sulA encoding SOS cell division inhibitor was induced by up to 
3.65 log2-fold from 30 min of the shift onward. This was 
accompanied by up-regulation of gene encoding the cell 
filamentation protein (fic) from 80 min onward (1.74 log2-fold on 
average). The products of both genes have previously been 
demonstrated to play an important role in cell filamentation 
[90,91]. Additionally, the apparent increase in suLA expression 
suggests that DNA damage might occur as a result of osmotic 
stress. This is in line with the previous findings of Fitt et al. [92] 
and Ribeiro et al. [93]. 

It is noteworthy that E. coli exposed to an hyperosmotic shift did 
not show a strong alteration in expression of genes and proteins 
involved in the biosynthesis of alarmone (p)ppGpp that instigates 
the stringent response [94] (Table S2 in File SI). This was despite 
that the proteomic data revealing a significant increase in 
abundance of SpoT (2.08 log2-fold) only at 310 min of the shift. 
However, previous studies have reported that the mechanism of 

References 

1. Cagnry C, Crowley H, Duffy G, Sheridan JJ, Brien SO, et al. (2004) Prevalence 
and numbers of Escfimchia coli 0157:H7 in minced beef and beefburgers from 
butcher shops and supermarkets in the Republic of Ireland. Food Microbiol 21; 

203-212. 

2. McClurc P (2000) The impact of £. coli 0157 on the food industry. 
WorldJ Microbiol Biotechnol 16: 749-755. 

3. Elder RO, Keen JE, Siragusa GR, Barkocy-Gallagher GA, Koohmaraie M, 
et al. (2000) Correlation of enterohemorrhagic Eschericfm cokOl^l prevalence in 
feces, hides, and carcasses of beef cattle during processing. Proc Nad Acad 
Sci U S A 97: 2999-3003. 

4. Weber A, Kogl SA, Jung K (2006) Time-dependent proteome alterations under 
osmotic stress during aerobic and anaerobic growth in Escherichia coli. J Bacteriol 
188: 716,5-7175. 

5. Csonka LN (1989) Physiological and genetic responses of bacteria to osmotic 
stress. Microbiol Rev 53: 121-147. 

6. Chung HJ, Bang \V, Drake MA (2006) Stiess response oi Escherichia coli, Compr 
Rev Food Sci Food Saf 5: 52-64. 

7. O'Byrne CP, Booth IR (2002) Osmoregulation and its importance to foodbome 
microorganisms. Int J Food Microbiol 74: 203-216. 

8. Gunasekera TS, Csonka LN, Paliy O (2008) Genome-wide tianscriptional 
responses of Escherichia coli K-12 to continuous osmotic and heat stresses. 
J Bacteriol 190: 3712-3720. 



(p)ppGpp action is typically characterized by a major change in 
global expression with the negatively controlled elements involved 
in cell growth (i.e. translation processes) and the positively 
controlled elements involved in amino acid biosynthesis, tran- 
scriptional factors, or alternative sigma factors (i.e. RpoE, RpoH 
and RpoS) [94,95] . This indeed corresponds well with the present 
observation on the response of E. coli to hyperosmotic stress as 
described above, suggesting that the stringent response might be 
induced under such stress. 

Conclusion 

The present study provides a deeper insight into the growth 
kinetics and potential mechanisms enabling survival and growth of 
E. coli 0157:H7 Sakai when subjected to a sudden downshift in 
water activity. Such knowledge wUl be used to develop more 
targeted approaches for the food industry to eliminate or control 
this pathogen. 

Supporting information 

File SI Table SI Summary on % false-positive discovery rate 
calculated for each MudPIT run. Table S2 Comparison of 
differentially expressed genes and proteins in Escherichia coli 
0157:H7 Sakai during exposure to hyperosmotic shift. Table S3 
T-profiler analysis of transcriptomic and proteomic data based on 
JCVI CMR functional categories. Table S4 T-profiler analysis of 
transcriptomic and proteomic data based on Ecocyc and KEGG 
functional categories. Table S5 T-profiler analysis of transcrip- 
tomic and proteomic data based on major regulons. 
(XLSX) 

Acl(nowledgments 

Wc would like to thank Ian Jcnson, Manager of the Food Salcty and 
Strategic Science program from Meal and Livestock Australia lor his 
support and guidance on the project. We also wish to thank Edwin Lowe 
and Richard Wilson from Central Science Lalioratory (University of 
Tasmania, Hobart, Australia) for expert technical support in 2D-LC/MS/ 
MS analysis. 

Author Contributions 

Conceived and designed the experiments: TR CK TK. Performed the 
experiments: CK TK. Analyzed the data: CK TK. Contributed reagents/ 
materials/analysis tools: CK. Wrote the paper: CK TK KG JPB TR. 



9. Kocharunchitt C, King T, Gobius K, Bowman JP, Ross T (2011) Integrated 
transcriptomic and proteomic analysis of the physiological response of Escherichia 
coli 0157:H7 Sakai to steady-state conditions of cold and water activity stress. 
Mol Cell Proteomirs. doi:lo'l074/mcp.Mll 1 1.009019. 

10. Hayashi T, Makino K, Ohnishi M, Kurokawa K, Ishii K, ct al. (2001) Complete 
Genome Sequence of enterohemorrphagic Eschedchia coli ()157:H7 and genomic 
comparison vrith a laboratory strain K-12. DNA Res 8: 1 1-22. 

11. Watanabe H, Wada A, Inagaki Y, Itoh K, Tamura K (1996) Outbreaks of 
enterohaemorrhagic Eschm£hia co/z 0157;H7 infection by two different genotype 
strains in Japan. Lancet 348: 831-832. 

12. Chiiife J, Resnik SL (1984) Unsaturated solutions of sodium chloride as 
reference sources of water activity at various temperatures. J Food Sci 49: 1486— 
1488. 

13. Mellefont LA, McMeekin TA, Ross T (2003) The cifcct of abrupt osmotic shifts 
on the lag phase duration of foodbome bacteria. IntJ Food Microbiol 83: 281— 
293. 

14. Robinson TP, Ocio MJ, Kaloti A, Maekey BM (1998) The effect of growth 
environment on the lag phase of Listeria tmmcytogenes. Int J Food Microbiol 44: 
83-92. 

15. Mellefont LA, Ross T (2003) The effect of abrupt shifts in temperature on the lag 
phase duration of Escherichia coli and Klebsiella oxytoca. Int J Food Microbiol 83: 
295-305. 



PLOS ONE I www.plosone.org 



11 



March 2014 | Volume 9 | Issue 3 | e90422 



E. coli Physiology under Water Activity Stress 



16. King T, Kocharunchitt C, Gobius K, Bowman J P, Ross T (2013) Global 
genome response of Escherichia coli 0157:H7 Sakai during dynamie ehangcs in 
growth kinetics induced by an abrupt temperature downshift. 

17. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, ct al. (2003) 
Exploration, normalization, and summaries of high density oHgonucleotide array 
probe level data. Biostatistics 4: 249^264. 

If5. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical 
and powerful approach to multiple testing. J R Stat Soc, Series B 57: 289—300. 

19. Rauch A, Bellew M, Eng J, Fitzgibbon M, Holzman T, et al. (2006) 
Gomputational proteomics analysis system (CPAS): an extensible, open-source 
analytic system for evaluating and publishing proteomic data and high 
throughput biological experiments. J Proteome Res 5: 112—121. 

20. Craig R, Bcavis RG (2004) TANDEM: matching proteins with tandem mass 
spectra. Bioinformatics 20: 1466-1467. 

21. Nesvizhskii Al, Keller A, Kolker E, Aebersold R (2003) A statistical model for 
identifying proteins by tandem mass speetrometr\'. Anal Chem 75: 4646—4658. 

22. Keller A, Nesvizhskii Al, Kolker E, Aebersold R (2002) Empirical statistical 
model to estimate the accuracy of peptide identifications made by MS/MS and 
database search. Anal Ghem 74: 5383-5392. 

23. Zybailov BE, Mosley AE, Sardiu ME, Coleman MK, Florcns E, et al. (2006) 
Statistical analysis of membrane proteome expression changes in Saccharomyces 
ceremsiae. J Proteome Res 5: 2339-2347. 

24. Pham TV, Piersma SR, Warmoes M, Jimenez CR (2010) On the beta-binomial 
model for analysis of spectral count data in label-free tandem mass spectrometry- 
based proteomics. Bioinformatics 26: 363—369. 

25. Boorsma A, Foat B(^, Vis D, Klis F, Bussemaker HJ (2005) T-profiler: scoring 
the activity of predefined groups of genes using gene expression data. Nucleic 
Acids Res'33: w592-w595. 

26. Price NE, Raivio TL (2009) Characterization of the Cpx regulon in Escherichia 
co/i" strain MC4 100. J Bacteriol 191: 1798-1815. 

27. Rhodius VA, Suh WC, Nonaka G, West J, Gross CA (2006) Conserved and 
v£U"iable fimctions of the stress response in related genomes. PLoS Biol 4: e2. 

28. Nonaka G, Blankschien M, Herman C, Gross CA, Rhodius VA (2006) Regulon 
and promoter analysis of the E. coli heat-shock factor, a"^^, reveals a multifaceted 
cellular response to heat stress, (jcnes Dev 20: 1776-1789. 

29. Patten GE, Kirchhof JVKi, Schertzberg MR, Morton RA, Schellhorn HE (2004) 
Microarray analysis of RpoS-mediated gene expression in Escherichia to/i K-12. 
Mol Genet Genomics 272: 580-591. 

30. Lacour S, Landini P (2004) c^-dependent gene expression at the onset of 
stationary phase in Escherichia coli: function of a -dependent genes and 
identification of their promoter sequences. J Bacteriol 186: 7186-7195. 

3 1 . Weber H, Polen T, HeuveHng J, Wendisch VF, Hengge R (2005) Genome-wide 
analysis of the general stress response network in Escherichia coli: CT^-dependent 
genes, promoters, and sigma factor selectivity. J Bacteriol 187: 1591—1603. 

32. Hagiwara D, Sugiura M, Oshima 'E, Mori H, Aiba H, et al. (2003) Genome- 
wide analyses re\ ealing a signaling network of the RcsG-YojN-RcsB phosphor- 
clay system in Escherichia coli.] Bacteriol 185: 5735—5746. 

33. de Hoon MJ, Imoto S, Nolan J, Miyano S (2004) Open source clustering 
software. Bioinformatics 20: 1453-1454. 

34. Saldanha AJ (2004) Java Treeview-extensible visuadization of microarray data. 
Bioinformatics 20: 3246-3248. 

35. Rustici G, Kolesnikov N, Brandizi M, Burdett T, Dylag M, et al. (2013) 
ArrayExpress update — trends in database growth and links to data analysis tools. 
Xueleir Acids Res 11: D987 D990. 

36. Vi/caino jA, Cotr R(;, Gsordas A, Dianes JA, Eabregat A, et al. (2013) The 
proteomics idrnlilications (PRIDE) database and associated tools: status in 2013. 
Nucleic Acids Res 41: D1063-D1069. 

37. Ross T, Ratkowsky DA, Mcllefont EA, McMcckin TA (2003) Modelling the 
effects of temperature, water activity, pH and lactic acid concentration on the 
growth rate of Escherichia coli. Int J Food Microbiol 82; 33-43. 

38. Mellefont LA, McMeekin TA, Ross T (2005) Viable count estimates of lag time 
responses for Salmonella typhimurium M48 subjected to abrupt osmotic shifts. 
Int J Food Microbiol 105: 399-410. 

39. Robinson TP, Aboaba GO, Kaloti A, Ocio MJ, Baranyi J, ct al. (2001) The 
effect of inoculum size on the lag phase of Listeria monocytogenes. Int J Food 
Microbiol 70: 163-173. 

40. Busta FF (1978) Introduction to injury and repair of microbial cells. Adv Appl 
Microbiol 23: 195-201. 

41. Zhou K, George SM, Metris A, Li PL, Baranyi J (2010) Lag phase of Salmonella 
under osmotic stress Appl Environ Microbiol 77: 1758-1762. 

42. Mackey BM, Derrick CM (1982) The effect of subletiial injury by heating, 
freezing, drying and gamma-radiation on the duration of the lag phase of 
Salmonella tpphimurium. J Appl Microbiol 53: 243-251. 

43. Chen G, Gharib TG, Huang CC, Taylor JM, Misek DE, et al. (2002) 
Discordant protein and mRNA expression in lung adenocarcinomas. Mol Cell 
Proteomics 1; 304-313. 

44. Ciriliin IJ, (iygi SP, Ideker T, Rist B, Eng J, et al. (2002) Complementary 
proliling of gene expression at the transcriptome and proteome levels in 
Saccharomyces cerevisiae. Mol Cell Proteomics 1: 323-333. 

45. Gygi SP, Rochon Y, Franza BR, Aebersold R (1999) Correlation between 
protein and mRNA abundance in yeast. Mol Cell Biol 19: 1720-1730. 

46. Rossouw D, van den Dool AH, Jacobson D, Bauer FF (2010) Comparative 
transcriptomic and proteomic proflEng of industrial wine yeast strains. Appl 
Environ Microbiol 76: 3911-3923. 



47. Hatzimanikatis V, Ghoe EH, Eee KH (1999) Proteomics: theoretical and 
experimental considerations. Biotechnol Prog 15: 312-318. 

48. Hatzimanikatis V, Lee KH (1999) Dynamical analysis of gene networks requires 
bodi mRNA and protein ejqjression information. Metab Eng 1: 275-281. 

49. Gyaneshwar P, P^diy O, McAuIiffe J, Jones A, Jordan MI, et al. (2005) Lessons 
fi-om Escherichia coli genes similarly regulated in response to nitrogen and sulfiir 
limitation. Proc Natl Acad Sci U S A 102: 3453-3458. 

50. Klumpp S, Zhang Z, Hwa T (2009) Growth rate-dependent global effects on 
gene expression in bacteria. Cell 139: 1366-1375. 

51. Weber A, Jung K (2002) Proliling early osmostress-dependent gene expression in 
Escherichia coli using DNA macroarrays. J Bacteriol 184: 5502—5507. 

52. Pin C, Rolfe MD, Munoz-Cuevas M, Hinton JC, Peck MW, et al. (2009) 
Network analysis of the transcriptional pattern of young and old cells of 
Escherichia coli during lag phase. BMC Syst Biol 3. 

53. Meury J (1994) Immediate and transient inhibition of the respiration of 
Escherichia coli under hyperosmotic shock. FEMS Microbiol Lett 121: 281—286. 

54. Aertsen A, Miehiels GW (2004) Stress and how bacteria cope with death and 
survival. Grit Rev Microbiol 30: 263-273. 

55. Zhou E, Eei X-H, Bochner BR, Wanner BE (2003) Phenotype microarray 
analysis of Escherichia coli K-12 mutants with deletions of all two-component 
systems. J Bacteriol 185: 4956-^972. 

56. Pratt LA, Hsing WH, Gibson KE, Silhavy TJ (1996) From acids to osm^. 
multiple factors influence synthesis of the OmpF and OmpC porins in Escherichia 
coli. Mol Microbiol 20: 911-917. 

57. Brenner M, Ames BN (1971) The histidine operon and its regulation. In: H.J 
Vogel, editor editors. Metabolic pathways. New York: Academic Press Inc. 

58. CasadesusJ, RothJR (1989) Absence of insertions among spontaneous mutants 
oi Salmonella typhimurium. Mol Gen Genet 216: 210—216. 

59. Frandscn N, D'Ari R (1993) Excess histidine enzymes cause AICAR- 
independent filamentation in Escherichia coli. Mol Gen Genet 240: 348-354. 

60. Ron EZ, Alajem S, Biran D, Grossman N (1990) Adaptation of Escherichia coli to 
elevated temperatures: the metA gene product is a heat shock protein. Antonie 
Van Leeuwenhoek 58: 169—174. 

61. Gur E, Biran D, Ga^it E, Ron EZ (2002) In vivo aggregation of a single enzyme 
limits growth of Escherichia coli at elevated temperatures. Mol Microbiol 46; 
1391-r397. 

62. Horinouchi 'E, Tamaoka K, Furusaw'a G, Gno N, Suzuki S, et al. (2010) 
Transcriptome analysis of parallel-evolved Escherichia coli strains under ethanol 
stress. BMG Genomics 11. 

63. Wood JM (1999) Osmosensing by bacteria: signals and membrane-based 
sensors. Microbiol Mol Biol Rev 63: 230-262. 

64. Jung IL, Oh TJ, Kim IG (2003) Abnormal growth of polyamine-deficient 
Escherichia coli mutant is partially caused by oxidative stress-induced damage. 
Arch Biochem Biophys 418: 125-132. 

65. Sleator RD, Hill C (2001) Bacterial osmoadaptation: the role of osmolytes in 
bacterial stress and virulence. FEMS Microbiol Rev 26: 49—71. 

66. Russell NJ, Evans Rl, ter Steeg PF, Hellemons J, Verheul A, et al. (1995) 
Membranes as a target lor stress adaptation. Int J Food Microbiol 28: 255-261. 

67. Canet S, Heyde M, Portalier R, Ealoi P (2003) Involvement of phospholipids in 
resistance and adaptation of Escherichia coli to acid conditions and to long-term 
survival. FEMS Microbiol Lett 225: 207-21 1. 

68. Chang YY, Cronan JE (1999) Membrane cyclopropane fatty acid content is a 
major factor in acid resistance of Escherichia coli. Mol Microbiol 33: 249-259. 

69. Shabala L, Ross T (2008) Cyclopropane fatty acids improve Escherichia coli 
survival in acidified minimal media by reducing membrane permeability to H 
and enhanced ability to extrude H^. Res Microbiol 159: 458—461. 

70. Brown JE, Ross V, McMeekin 'EA, Nichols PD (1997) Acid habituation of 
Escherichia coli and the potential role of cyclopropane latty acids in low pH 
tolerance. IntJ Food Microbiol 37: 163-173. 

71. Guillot A, Obis D, Mistou MY (2000) Fatty acid membrane composition and 
activation of glycine-betaine transport in Lactococcus lactis subjected to osmotic 
stress. IntJ Food Microbiol 55: 47-51. 

72. Vaara M, Nurminen M (1999) Outer membrane permeability barrier in 
Escherichia coli mutants that are defective in the late acyltransferases of lipid A 
biosynthesis. Antimicrob Agents Chemother 43: 1459—1462. 

73. Vuorio R, Vaara M (1992) 'Ehe lipid A biosynthesis mutation lpxA2 of Escherichia 
coli results in drastic antibiotic supersusceptibility. Antimicrob Agents Che- 
mother 36: 826-829. 

74. Meredith TC, Mamat U, Kaczynski Z, Einder B, Hoist O, et al. (2007) 
Modification of lipopolysaccharide with colanic acid (M-anti^en) repeats in 
Escherichia co£ J Biol Chem 282: 7790-7798. 

75. Majdalani N, Heck M, Stout V, Gottesman S (2005) Role of ResF in signaling to 
the Res phosphorelay pathway in Escherichai coli.] Bacteriol 187: 6770 6778. 

76. Ophir T, Gutnick DE (1994) A role for exopolysaecharides in the protection of 
mieroorgansims from desiccation. Appl Environ Microbiol 60: 740—745. 

77. Chen J, Lee SM, Mao Y (2004) Protective effect of cxopolysaccharide colanic 
acid of Escherichia coU 0157:H7 to osmotic and oxidative stress. IntJ Food 
Microbiol 93: 281-286. 

78. Mao YM, Doyle P, Chen J (2006) Role of colanic acid exopolysaccharide in the 
survival of enterohaemorrhagic Escherichia coli 0157:H7 in simulated gastroin- 
testinal fluids. Lett Appl Microbiol 42: 642-647. 

79. Danese PN, Pratt EA, Kolter R (2000) Exopolysaccharide production is rcc^uired 
for development of Escherichia coli K-12 biofilm architecture. J Bacteriol 182: 
3593-3596. 



PLOS ONE I www.plosone.org 



12 



March 2014 | Volume 9 | Issue 3 | e90422 



E. coli Physiology under Water Activity Stress 



80. Sl^di^■^ki DL), Ciottcsman S (1996) Osmotic shock induction of capsule synthesis 
in Esdienchia colt K-12. J Bacteriol 178: 1204-1206. 

81. Allen PM, Fisher D, Saunders JR, Hart CA (1987) The role of capsular 
polysacchjiride K2 lb of Klebsiella and of the structurally related colanic-acid 
polysaccharide of Escherichia coli in resistance to phagocytosis and serum killing. 
J Med Microbiol 24: 363-370. 

82. Kim Y-H, Han KY, Lee K, LeeJ (2005) Proteome response of Escherichia to/z fed- 
batch culture to temperature downshift. Appl Microbiol Biotechnol 68: 786- 
793. 

83. Wecker P, Klockow C, EUrott A, Quast G, Langhammcr P, ct al. (2009) 
Transcriptional response of" the model planctomyccte Rhodopirellula baltka SHl""" 
to changing environmental conditions. BMC (Genomics 10: 410—425. 

84. Franccz-G harlot A, Laugel B, Gemcrt AV, Dubarry N, Wiorowski F, et al. 
(2003) RcsCDB His-Asp phosphorelay system negatively regulates the JlhDC 
operon in Escherichia coli. Mol Microbiol 49: 823-832. 

85. Zheng M, Wang X, Templeton LJ, Smulski DR, LaRossa RA, et al. (200 1) DNA 
microarray-mediated transcriptional profiling of the Escherichia coli response to 
hydrogen peroxide. J Bacteriol 183; 4562-4570. 

86. Alba BM, Gross CA (2004) Regulation of the Escherichia coli a^-dependent 
envelope stress response. Mol Microbiol 52: 613-619. 



87. Meury J, Kohiyama M (1991) Role of heat shock protein DnaK in osmotic 
adaptation of Escherichia coli.] Bacteriol 173: 4404—4410. 

88. Thomas JG, Baneyx F (2000) ClpB and HtpG facilitate de novo protein folding in 
stressed Escherichia coli cells. Mol Microbiol 36: 1360-1370. 

89. Bianchi AA, Baneyx F (1999) Hyperosmotic shock induces the and <T^ stress 
regulons of Escherichia coli. Mol Microbiol 34: 1029-1038. 

90. Bi E, Lutkenhaus J (1993) Gell division inhibitors SulA and MinGD prevent 
formation of the FtsZ ring. J Bacteriol 175: 1 1 18-1 125. 

91. Kawamukai M, Matsuda H, I'ujii W, Xishida T, Izumoto Y, et al. (1988) 
Gloning of the flc-\ gene in\'ol\'ed in cell filamcntation induced by cyclic AJMP 
and construction of a delta Jic Escherichia coli strain. J Bacteriol 170: 3864—3869. 

92. Fitt PS, Sharma N, BaU D (1992) Induction of error-free DNA repair in 
Escherichia coli by nonmutagenic stress. Curr Microbiol 24: 251-256. 

93. Ribeiro GF, Gorte-Real M, Johansson B (2006) Characterization of DNA 
damage in yeast apoptosis induced by hydrogen peroxide, acetic acid, and 
h^'pcrosmoite shock. Mol Biol Gell 17: 4584-4591. 

94. Traxler MF, Summers SM, Nguyen H-T, Zacharia VM, Hightower (iA, ct al. 
(2008) The global, ppGpp-mediatcd stringent response to amino acid starvation 
in Escherichia coli. Mol Microbiol 68: 1 128-1 148. 

95. Durfee T, Hansen A-M, Zhi H, Blattner FR, Jin DJ (2008) Transcription 
profiling of the stringent response in Escherichia coli.] Bacteriol 190: 1084—1096. 



PLOS ONE I www.plosone.org 



13 



March 2014 | Volume 9 | Issue 3 | e90422 



